Electromagnetically Induced Transparency and Absorption in Metamaterials: 
The Radiating Two-Oscillator Model and Its Experimental Confirmation 
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Several classical analogues of electromagnetically induced transparency in metamaterials have 
been demonstrated. A simple two-resonator model can describe their absorption spectrum qualita- 
tively, but fails to provide information about the scattering properties — e.g., transmission and group 
delay. Here we develop an alternative model that rigorously includes the coupling of the radiative 
resonator to the external electromagnetic fields. This radiating two-oscillator model can describe 
both the absorption spectrum and the scattering parameters quantitatively. The model also predicts 
metamaterials with a narrow spectral feature in the absorption larger than the background absorp- 
tion of the radiative element. This classical analogue of electromagnetically induced absorption is 
shown to occur when both the dissipative loss of the radiative resonator and the coupling strength 
are small. These predictions are subsequently demonstrated in experiments. 
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Electromagnetically induced transparency (EIT) is an 
effect that renders an otherwise opaque medium trans- 
parent in a narrow transmission window with low absorp- 
tion and steep dispersion [T] . It has attracted quite some 
interest, because of its promise for a low-loss slow- light 
medium. EIT was first demonstrated in certain three- 
level atomic systems like alkali vapors, where destructive 
interference between two radiative transitions creates a 
dark dressed superposition state with no electric dipole 
moment [2H5]. Quantum- mechanical EIT allows for the 
slowdown of light to a group velocity of about 17 m/s [6 
and even for the storage of light [7H1], but it requires 
complicated experimental handling because of the rather 
short coherence times of the superposition state. 

However, it was realized soon that the characteris- 
tic features — simultaneously low absorption and steep 
dispersion — can also be realized in classical systems such 
as coupled mechanical or electrical resonators [TU] or even 
coupled acoustic resonators [IT] . This has led to the 
demonstration of many classical analogues of EIT, e.g., 
in electromagnetic metamaterials [12H23) and optical mi- 
croresonators [23H50]. A simple model to describe these 
systems is a set of two coupled harmonic oscillators, 
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The radiative resonator with resonance frequency w r and 
damping factor 7 r is described by the excitation p(t) and 
is driven by the external force f(t). The dark resonator 
with resonance frequency Wd and damping factor 7d is 
described by the excitation q(t). Both resonators are 
linearly coupled with coupling strength k. The individ- 
ual oscillators can be mechanical, molecular, or subwave- 
length electromagnetic elements. The excitations would 
then represent the corresponding physical quantities such 
as the displacement from the rest position (mechanical), 



or the microscopic electric or magnetic dipole moment. 
This model reflects the essential ingredients of EIT: two 
coupled resonances that are are asymmetrically driven 
by the external force. Eqs. ([I]) and ^ can be solved in 
the frequency domain by assuming a solution of the form 
p(t) = p{u>) exp(— iuit) and q{t) = q(u>) exp(— hut): 
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where D T (u>) = 1 — (w/w r ) 2 — vy T (u)/uj T ) and = 
1 — (w/ajd) 2 — i7d(w/wd). The dissipated power per unit 
cell, which can be obtained from Q = ^(7 r |p(a;)| 2 + 
7d|(?(cLi)| 2 ), has a Lorentzian shape with a sharp incision 
at the resonance frequency if uj t s» U) d , 7d -C 7 r , and 
7d7r < k 2 < 1. 

Even though the two-resonator model can qualitatively 
describe the absorption of classical EIT analogues, it fails 
to model scattering parameters of metamaterials exhibit- 
ing a classical EIT response. This is especially trouble- 
some since it makes it impossible to determine the group 
delay, quite an essential parameter for slow-light media. 

In this Letter, we will develop a slightly more complex 
model for classical EIT media. In contrast to the micro- 
scopic two-oscillator model, which does not contain infor- 
mation about the actual coupling to the external world, 
we introduce our new radiating two-oscillator model that 
rigorously describes both the microscopic and the macro- 
scopic response in terms of the radiated field (i.e., the 
incident, reflected, and transmitted waves). We start its 
derivation by recognizing that most of the EIT meta- 
materials fabricated to date are essentially single-layer 
structures rather than bulk media. Hence, their effective 
response can better be described by an electric current 
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sheet with surface conductivity a sc (we restrict the dis- 
cussion here to metamaterials with an electric dipole re- 
sponse, but the same analysis can be applied to metama- 
terials with a magnetic dipole response). The scattering 
parameters of an electric current sheet are [31] 



the absorbance and the group delay: 



R 



(4) 
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where £ is the wave impedance of the external waves. 
Equations Q serve as the world model; i.e., they de- 
scribe the interaction of the medium with the external 
electromagnetic field. The microscopic behavior of the 
EIT medium can still be described by the two-resonator 
model as given by Eqs. and (|2|). 

In order to complete the radiating two-oscillator 
model, we have to find a connection between the external 
behavior of the system (the surface field E s and the sur- 
face conductivity o~ se ) and the microscopic behavior (the 
excitations p and q and the driving force /). The macro- 
scopic surface field E s is the spatially averaged electric 
field on the current sheet — it is related to the incident 
field by E s = T x Ei n and can be calculated in this way 
both for experiments and simulations. First, we observe 
that each of the constituent meta-atoms contributes a 
dipole moment p to the metamaterial and, if there are n s 
atoms per unit of surface area, the average polarization 
current thus equals 



O's(i)) = n s p{t) o (jsM) = -ium s p(w). 



(5) 



The dark resonator does not contribute to the surface 
current since it has no dipole moment commensurate 
with the external field. Secondly, we need to find a con- 
nection between the surface field E s , which drives the 
dipole oscillation in the world model, and the driving 
force / in the microscopic model; i.e., we seek the pro- 
portionality constant C in f(t) — CE s (t). (Note that 
the surface field E s is different from the incident field be- 
cause of the scattering from the meta-atoms.) This can 
be done by recalling that for our linear meta-atom the 
average dipole moment must be proportional to the elec- 
tric field at the surface: n s p(uj) = e x so (w)i?s(w), where 
Xso is the surface susceptibility. In the static limit, this 
yields 

e xiT tic) E s (0) = n s p(Q) =n s (l- n 2 )' 1 /(0) « nj(0), 

(6) 

where we used the fact that K -C 1 under EIT condi- 
tions in the last approximation. Using Eqs. ^ and ([6]), 
we can now determine the surface conductivity from the 
constitutive equation (j s (oj)) = a se E s (u}): 
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where (3 = £oXsc tatIC ''- Once we have determined the sur- 
face conductivity, we can calculate the scattering param- 
eters from Eqs. Q and other derived quantities, such as 
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Note that, from the perspective of the microscopic de- 
scription in terms of two coupled oscillators used in pre- 
vious literature, our model introduces a radiation damp- 
ing term in the bright oscillator as well as an excitation- 
dependent external driving force, both of which originate 
in the scattered field of the bright resonators responsible 
for the reflectance (and nonunity transmittance) of the 
macroscopic sample. However, even with those correc- 
tions to the two-resonator model, we still need the full 
radiating two-oscillator model as described in this Letter 
to calculate the transmittance and group delay — the two 
most important characteristics of an EIT system. 

The radiating two-oscillator model allows us to under- 
stand the response of EIT metamaterials as a function 
of their microscopic parameters. In Fig. [lja), we plot 
the surface conductivity, the absorption, the transmis- 
sion amplitude and phase, and the group delay for a set 
of metamaterials with different dissipativc damping in 
the radiative resonator. In the top row (high damping), 
we recognize the typical features of EIT. The conduc- 
tivity has a Lorentzian envelope with a sharp incision, 
resulting in a frequency window with large transmission 
and small absorption. At the same time, there is large 
normal dispersion in the transmission phase, which leads 
to a significantly enhanced group delay. The reduced 
response can be understood from the destructive inter- 
ference of the excitation p due to the external field and 
due to the coupling from the dark resonator (see \p(<u))\ 
at the resonance frequency in Fig. SI [3"2"]). 

An interesting phenomenon occurs if we decrease the 
dissipative loss factor j T of the radiative resonator [sec- 
ond row in Fig. [lja)]. There is still a frequency window 
with high transmission, but the incision in the absorp- 
tion spectrum becomes smaller and finally disappears. 
This does not mean, however, that the dark resonance 
has disappeared, as we can clearly see from the strong 
dispersion and associated enhanced group delay. Rather, 
the background absorption of the radiative resonance is 
decreased, but the radiative resonance is still sufficiently 
broadened by the radiation damping, while the excita- 
tion of the dark resonance is barely changed. At a cer- 
tain point, the absorption reduction in the radiative res- 
onator is exactly cancelled by the absorption in the dark 
resonator. When we further decrease the dissipative loss 
of the radiative resonator, the absorption spectrum turns 
into a very dim background with a narrow peak at the 
resonance frequency oj = lu^. This phenomenon can be 
seen as a classical analogue of electromagnetically induced 
absorption (EIA) |33j and we believe it may be interest- 
ing for applications in spectroscopy and sensing since the 
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FIG. 1: (Color online) Spectra of the surface conductivity, absorption, transmission amplitude and phase, and group delay of 
EIT or EIA metamaterials as described by the radiating two-oscillator model, (a) As a function of the dissipative damping 
factor of the radiative resonator. (74 = 0.0015, k = 0.06, and /3 = 1.) (b) As a function of the dissipative damping factor of 
the dark resonator. (7 r = 0.1, K = 0.06, and ft — 1.) (c) As a function of the coupling strength. (7d = 0.0015, jr = 0.1, and 
= 1.) 



width of the peak is reduced by the lack of radiation 
damping in the dark resonator and the additional nar- 
rowing due to the coupling. Note that the EIA effect 
could only be described by the radiating two-oscillator 
model, since the bare two-resonator model lacks radia- 
tive broadening of the radiative resonator. 

The transition between EIT and EIA can also be ob- 
served when we increase the dissipative loss of the dark 
resonator [see Fig. [TJb)]. When 7^ is increased, the ra- 
diative resonance remains unaltered, but the absorption 
in the center of the transparency window goes up. Even- 



tually, the loss in the dark resonator overcomes the loss 
reduction due to the destructive interference in the ra- 
diative resonator. Note, however, that too large a value 
of 7d destroys the EIT or EIA phenomenon as shown in 
Fig. S2 [35]. Finally, EIA can also be achieved by de- 
creasing the coupling strength k, as in Fig.Qc). Weaker 
coupling creates a narrower transparency window with 
larger excitation q in the dark resonator. This in turn 
increases the absorption at the resonance frequency, re- 
sulting in EIA when the dissipation in the dark resonance 
overcomes the loss reduction in the radiative resonance. 
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FIG. 2: (Color online) Absorption spectra of the EIT or EIA metamaterials for different coupling strengths, (a) Predicted by 
the radiating two-oscillator model with parameters lj^ — 2n x 10.58 GHz, 7d = 0.0035, uu T — 2n x 10.12 GHz, 7 r = 0.01, and 
P = 0.33/(cj r r;o). (b) Measured for microwave metamaterials in the X band. 
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FIG. 3: (Color online) Absorption spectra of the EIT or EIA metamaterials for different damping of the radiative resonator, 
(a) Predicted by the radiating two-oscillator model with parameters uj r — 2n x 10.11 GHz, 7d = 0.0032, k = 0.027, and 
P = 0.38/(td r 77o). (b) Measured for microwave metamaterials in the X band. The inset shows the linear relationship between 
the damping constants in the model and the lumped resistances in the experiment. The error bars represent the 5% tolerance 
of the resistors. 



Again, EIT or EIA is destroyed when n 2 < j^Ji- (see 
Fig. S3 [32]). 

Subsequently, we have confirmed the predictions of the 
radiating two-oscillator model in microwave metamateri- 
als consisting of a copper cut wire as the radiative res- 
onator and two copper closed-ring resonators as the dark 
resonator on a Rogers substrate [223 ( see inset in Fig. [2]) . 
The closed rings provide a dark resonance by using the 
antisymmetric hybridization of the electric dipole modes 
of the split-ring resonators — this hybridization is actu- 
ally an electric quadrupole resonance with zero overlap 
with the fundamental waveguide mode (see the Supple- 
mental Material for field plots demonstrating the dark 
and bright modes [32 ). The measurements were per- 
formed in a WR-90 waveguide and the scattering param- 
eters were measured using a vector network analyzer (HP 
E8364) and calibrated using a transmission-reflection- 
match (TRM) method. In Fig. [2jb), we have plotted 
the experimental absorption spectrum as a function of 
the coupling strength (the coupling strength was varied 



by moving the cut wire over a distance d horizontally 
between the closed rings). We see that the absorption 
spectrum of the model [Fig. [2ja)] is in excellent agree- 
ment with the experiments [Fig. |2jb)]. Figure S5 shows 
that also the transmission and reflection coefficients are 
in excellent agreement [32 . The values of k from the 
model are perfectly proportional to the offset of the wire 
from the center (d), further confirming that the model 
can quantitatively describe EIT or EIA metamaterials. 
In a second set of experiments, we have altered the dissi- 
pative loss of the radiative resonator (by soldering SMD 
resistors with different resistance values R into the cut 
wire). We again see that the experiments [Fig. [3^b)] 
convincingly reproduce the transition from EIT to EIA, 
and are in very good agreement with the radiating two- 
oscillator model [Fig. [3](a)]. The transmission and re- 
flection coefficients from the model are also in excellent 
agreement with the experiments (see Fig. S6 [32J). The 
matching damping constants are plotted as a function of 
the resistor values in the inset of Fig. [3j revealing a lin- 
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ear relationship, though with an offset this time. The 
offset is because part of the dissipation in the radiative 
resonator happens in the capacitor due to relaxation loss 
in the substrate. The equivalent resistance of the relax- 
ation loss is estimated to be 3.4 fL The experimental 
group delay curves are provided in Fig. S7 [52] . 

In this Letter, we have focused on metamaterials with 
subwavelength constituents. Our model can not only 
quantitatively describe EIT metamaterials, but it also re- 
veals a classical analogue of EIA — which is characterized 
by a sharp absorption peak on a shallow background — 
when the radiative resonator has small dissipative loss, 
but is still sufficiently broadened by radiation damping. 
There have recently been two papers in which a phe- 
nomenon similar to EIA is discovered when the two res- 
onators are coupled to the external wave with different 
phase |35| or when a retardation-induced phase shift oc- 
curs in the coupling mechanism [34]. We believe that 
our model can also describe these experiments with small 
changes (e.g., complex k). Nevertheless, in truly homog- 
enizable metamaterials, such phase differences are not 
possible and we must revert to radiation-broadened res- 
onators to achieve EIA. 
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